Enumeration of spanning trees in a pseudofractal scale-free web 
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Abstract. - Spanning trees are an important quantity characterizing the reliability of a network, 
however, explicitly determining the number of spanning trees in networks is a theoretical challenge. 
In this paper, we study the number of spanning trees in a small-world scale-free network and obtain 
the exact expressions. We find that the entropy of spanning trees in the studied network is less 
than 1, which is in sharp contrast to previous result for the regular lattice with the same average 
degree, the entropy of which is higher than 1. Thus, the number of spanning trees in the scale-free 
network is much less than that of the corresponding regular lattice. We present that this difference 
lies in disparate structure of the two networks. Since scale-free networks are more robust than 
regular networks under random attack, our result can lead to the counterintuitive conclusion that 
a network with more spanning trees may be relatively unreliable. 



Introduction. — The enumeration of spanning trees 
in networks (graphs) is a fundamental issue in mathemat- 
ics PUS] , physics [4]|5] , and other discipline [6] . A spanning 
tree of any connected network is defined as a minimal set 
of edges that connect every node. The problem of span- 
ning trees is relevant to various aspects of networks, such 
as reliability 0H], optimal synchronization [9], standard 
random walks [TUJ, and loop-erased random walks 
In particular, the number of spanning trees corresponds 
to the partition function of the g-state Potts model [13] 
in the limit of q approaching zero, which in turn closely 
relates to the sandpile model [13]. 

Because of the diverse applications in a number of 
fields P3] , a lot of efforts have been devoted to the study of 
spanning trees. For example, the exact number of span- 
ning trees in regular lattices [3HH] and Sierpinski gas- 
kets [5] has been explicitly determined in previous studies. 
However, regular lattices and fractals cannot well mimic 
the real-life networks, which have been recently found 
to synchronously exhibit two striking properties: scale- 
free behavior [16] and small-world effects [17] that has a 
strong impact on the enumeration problems on networks. 
For example, previous work on counting subgraphs, such 
as cliques [18], loops and Hamiltonian cycles [191 . has 
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shown that scale-free degree distribution implies a very 
non-trivial structure of subgraphs. However, so far in- 
vestigation on the number of spanning trees in scale-free 
small-world networks has been still missing. In view of 
the distinct structure, as compared to regular lattices, it 
is of great interest to examine spanning trees in scale-free 
small- world networks. 

In this paper, we intend to fill this gap by providing 
a first analytical research of spanning trees in a small- 
world network with inhomogeneous connectivity. In order 
to exactly obtain the number of spanning trees, by using 
a renormalization group method [50] , we consider a deter- 
ministically growing scale-free network with small-world 
effect. We find that the entropy of its spanning trees is 
smaller than 1, which is a striking result that is qualita- 
tively different from that of two dimensional regular lat- 
tices with identical average degree, in which the entropy 
is higher than 1. Thus, the number of its spanning trees is 
much lower than that of its corresponding regular lattice. 
We show that this difference can be accounted for by the 
heterogeneous structure of scale- free networks. Since the 
network under study is much robust to random deletion 
of edges, as opposed to regular lattice, our result suggests 
that networks with more spanning trees are not always 
more stable to random breakdown of edges, compared with 
those networks with less spanning trees. 
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Fig. 1: The first three generations of the iterative scale- free 
network. 



Pseudofractal scale-free web. — The studied scale- 
free network [5TJ[32], denoted by G n after n (n > 0) gen- 
erations, is constructed as follows: For n = 0, Go is a 
triangle. For n > 1, G n is obtained from G„_i: every 
existing edge in G„_i introduces a new node connected to 
both ends of the edge. Figure [T] illustrates the construction 
process for the first three generations. The network ex- 
hibits some typical properties of real networks. Its degree 
distribution P(k) obeys a power law P(k) ~ fci+lns/h^ 
the average distance scales logarithmically with network 
order (number of nodes) [23] , and the clustering coefficient 
is | . Alternatively, the network can be also created in an- 
other method [23|I24J . Given the generation n, G n +i may 
be obtained by joining at the hubs (the most connected 
nodes) three copies of G n , see Fig. [2] According to the 
latter construction algorithm, we can easily compute the 
network order of G n is V n = - — ^ n *-* n , there are 
three hubs denoted by A n , B n , and G„, respectively. 

Number of spanning trees. — After introducing the 
network construction and its properties, next we will study 
both numerically and analytically spanning trees in this 
scale-free network. 

Numerical solution. According to the well-known re- 
sult |25| , we can obtain numerically but exactly the num- 
ber of spanning trees, Nst(ti), by computing the non-zero 
eigenvalues of the Laplacian matrix associated with G n as 
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where Xi(n) (i = 1, 2, . . . , V n — 1) are the V n — 1 nonzero 
eigenvalues of the Laplacian matrix for G n . For a network, 
the non-diagonal element iy (i j) of its Laplacian matrix 
is -1 (or 0) if nodes i and j are (or not) directly connected, 
while the diagonal entry la equals the degree of node i. 
Using Eq. ([T]), we can calculate directly the number of 
spanning trees Nst{it>) of G„ (sec Fig. [3]). From Fig. [31 we 
can see that NsT{n) approximately grows exponentially in 
V n . This allows to define the entropy of spanning trees for 
G n as the limiting value [THS] 
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Fig. 2: (Color online) Second construction method of the net- 
work. G n +i can be obtained by joining three copies of G n 
denoted as g! 77 ' (r/ — 1, 2, 3), the three hubs of which are rep- 
resented by A^ , B^ , and Cffl . In the merging process, hubs 



A^ (resp. C y n>, A)t>) and B^' (resp. B)t\ C y n>) are identified 
as a hub node A n+ i (resp. B n+1 , C n +i) in G n+ i. 
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which is a finite number and a very interesting quantity 
characterizing the network structure. 

It should be mentioned that although the expression 
of Eq. ([1]) seems compact, the computation of eigenval- 
ues of a matrix of order V n x V n makes heavy demands 
on time and computational resources for large networks. 
Thus, one can count the number of spanning trees by di- 
rectly calculating the eigenvalues only for the first several 
iterations, which is not acceptable for large graphs. Par- 
ticularly, by using the eigenvalue method it is difficult and 
even impossible to obtain the entropy Ec n - It is thus of 
significant practical importance to develop a computation- 
ally cheaper method for enumerating spanning trees that 
is devoid of calculating eigenvalues. Fortunately, the itera- 
tive network construction permits to calculate recursively 
Nst(ji) and Ec n to obtain exact solutions. 

Closed-form formula. To get around the difficulties of 
the eigenvalue method, we use an analytic technique based 
on a decimation procedure [5D|. For simplicity, we use t n 
to express A^st(^)- Moreover, let a n denote the number 
of spanning subgraphs of G n consisting of two trees such 
that the hub node A n belongs to one tree and the two 
other hubs (B n and C n ) are in the other tree. Analo- 
gously, we can define quantities b n and c„, see Fig.[3J By 
symmetry, we have a n = b n = c„. Thus, in the following 
computation, wc will replace b n and c n by a n . 

Considering the self-similar network structure, the fol- 
lowing fundamental relations can be established: 

tn+l = (tn) 2 {a>n + c n + a„ + b n + c n + b n ) = 6a n (t n ) 2 (3) 



and 



a n +i = t n [(c n ) 2 + a n b n + b n c n + a n c n ] = At n {a n ) 2 . (4) 

Equation ([3J can be explained as follows. Since G„+i is 
obtained via merging three G n by identifying three couples 
of hub nodes, to get the number of spanning trees t n+ \ for 
G„+i, one of the copies of G n must be spanned by two 
trees. There are six possibilities as shown in Fig. [5[ from 
which it is easy to derive Eq. ([3]). Analogously, Eq. ([4} 
can be understood based on Fig. [6[ 
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Fig. 3: Logarithm of the number of spanning trees NsT(n) in 
network G n as a function of network order V n on a log-log 
scale. In the figure, L n = lnA r sx(n-); the filled symbols are 
the numerical results obtained from Eq. {l}, while the empty 
symbols correspond to the exact values from Eq. (|lip . both of 
which completely agree with each other. 



To obtain t n , we define an intermediary variable h T 
I s - that obeys the following recursive relation 
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With the initial condition to — 3 an( l a o = 1, we have 
ho = 3. Hence, Eq. §5§ is solved to yield 
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Plugging this expression into Eq. ([3]) leads to 
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Considering the initial value to = 3, we can solve Eq. (JSJ 
to obtain the explicit solution 

iVsx(n) = t n = 2(3" + 1 -2«-3)/4 3 (3"+ 1 +2 n+ l)/4 _ (g) 

Analogously, we can derive the exact formula for dji RSI 

On = 2 (3" +1 + 2 »-3)/ 4 3( 3 " +1 - 2 "- 3 )/4 . (10) 

It not difficult to represent Nst(ii) as a function of the 
network order V n , with the aim to obtain the relation be- 
tween the two quantities. Recalling V„ = - — j"^, we have 
3" +1 = 2V n - 3 and n + 1 = m(2V n - 3)/ In 3. These 
relations enable one to write Nst(ji) in terms of V n as 

N ST (n) = 2[ y «- ln (2^-3)/ln3-2]/2 3 H/„+ln(2V„-3)/ln3-2]/2 _ / 

We have confirmed the closed-form expressions for 
NsTin) against direct computation from Eq. (JTJ. In the 
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Fig. 4: Illustrative definition for the spanning subgraphs of 
G n . The two hub nodes connected by a solid line are in one 
tree, and the two hub nodes linked by a dotted line belong to 
different trees. 



full range of < n < 8, they are perfectly consistent with 
each other, which shows that the analytical formulas pro- 
vided by Eqs. © and pT|) are right. Figure [3] shows the 
comparison between the numerical and analytical results. 

Equation unveils the explicit dependence relation 
of NsT(n) on the network order V n . Inserting Eq. (|TT1) 
into Eq. ((2J, it is easy to obtain the entropy of spanning 
trees for G n given by 



lniV ST (n) 1, , . 

E Gn = Jim = -(In2 + ln3) ~ 0.89588 . (12) 
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This obtained asymptotic value is the smallest entropy 
(lower than 1) that has not been reported earlier for other 
networks with an average degree of 4. For example, the en- 
tropy for spanning trees in the square lattice is 1.16624 
a value larger than 1. Thus, the number of spanning trees 
in G n is much less than that in the square lattice with the 
same average degree of nodes. 

From the result obtained above, we can conclude that 
the pscudofractal scale-free network has less spanning 
trees than the regular lattice with the same average de- 
gree. The difference can be attributed to the structural 
characteristics of the two classes of networks. In scale- 
free networks, nodes have a heterogeneous connectivity, 
which leads to an inhomogeneous distribution of Lapla- 
cian spectra [21,26,27 . On the contrary, in regular lat- 
tices, since all nodes have approximately the same degree, 
their Laplacian spectra have a homogenous distribution. 
Thus, for two given scale-free and regular networks with 
the same order and average node degree, the sum of the 
eigenvalues of their Laplacian matrices are the same, but 
the product of non-zero Laplacian spectra of the scale- 
free network is smaller than its counterpart of the regu- 
lar network, because of the different distributions of the 
Laplacian spectra resulting from their distinct connectiv- 
ity distribution. Hence, the heterogeneous structure is 
responsible for the difference of number of spanning trees 
in scale-free networks and regular lattices. It should be 
stressed that although we only study a specific determin- 
istic scale-free network, we expect to find a qualitatively 
similar result about spanning trees in real-world scale- free 
Ijnetworks, since they have similar structural characteristics 
as that discussed above. 

As an important invariant of a network, the number of 
spanning trees is a relevant measure of the reliability of the 
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Fig. 5: (Color online) Illustration for the recursion expression 
for the number of spanning trees t n +i in G n +i- The two nodes 
at both ends of a solid line are in one tree, while the two nodes 
at both ends of a dotted line are in separate trees. 

network. Intuitively, among all connected graphs with the 
same numbers of nodes and edges, networks having more 
spanning trees are more resilient (reliable) to the random 
removal of edges, compared with those with less span- 
ning trees. That is to say, the former has a larger thresh- 
old of bond percolation than that of the latter. However, 
recent work [28H30] have shown that inhomogeneous net- 
works, such as scale- free networks, are impressively robust 
than homogeneous networks (e.g., exponential networks 
and regular networks) with respect to random deletion of 
edges. Thus, combining with our above result, we can 
reach the following counterintuitive conclusion that net- 
works (e.g., scale- free networks) with less spanning trees 
do not mean more vulnerable to random breakdown of 
links than those (e.g., regular lattices) with more span- 
ning trees. 

Conclusions. — In summary, diverse real-life net- 
works possess power-law degree distribution and small- 
world effect. In this paper, we have studied and enumer- 
ated explicitly the number of spanning trees in a scale- free 
network with small-world behavior. The exact solution 
was obtained on the basis of some precise recursion rela- 
tions derived from the iterative construction of the net- 
work addressed. It was demonstrated that scale-free net- 
work has much less spanning trees compared to the regu- 
lar lattice with the same number of nodes and edges. It 
was shown that this difference is rooted in the inherent 
architecture of the two types of networks. Although it is 
generally thought that increasing the number of spanning 
trees over all networks with identical number of nodes and 
edges can lead to a less fragile network, our results strik- 
ingly indicate otherwise. Our work may be helpful for 
designing and improving the reliability of networks. 
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Fig. 6: (Color online) Illustration for the recursive expression 
for the number of spanning subgraphs a n +i corresponding to 
network G n +i- The two nodes at both ends of a solid line 
(dotted line) are in one tree (two trees). 
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